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Abstract 

> 

Static, spherically symmetric solutions of the Yang-Mills-Dilaton theory are studied. It is shown 
that these solutions fall into three different classes. The generic solutions are singular. Besides 
, there is a discrete set of globally regular solutions further distinguished by the number of nodes 

of their Yang-Mills potential. The third class consists of oscillating solutions playing the role of 
limits of regular solutions, when the number of nodes tends to infinity. We show that all three 
sets of solutions are non-empty. Furthermore we give asymptotic formulae for the parameters 
q-i of regular solutions and confront them with numerical results. 

> : 

X ' 1 Introduction 



The dilaton may be considered as a kind of scalar graviton sharing with it a universal coupling 
to matter. From this point of view it may be not too surprising that the static, spherically 
symmetric solutions of the Yang-Mills-Dilaton (YMD) theory share many properties with their 
Einstein- Yang- Mills (EYM) relatives. In fact, numerical studies ^ |5] have revealed a great 
similarity between a family of 'gravitational sphalerons' - the Bartnik-McKinnon (BK) solutions 
|3j - and a corresponding family of dilaton solutions. Also an existence proof of these solutions 
running exactly along the lines of the one for the BK solutions jl] could be given [3]. 

Introducing a 'stringy' radial variable rescaled with a dilaton factor the field equations of 
the YMD theory resemble very much those of the EYM theory. At first sight the only difference 
is the larger number of gravitational degrees of freedom, which is however upset by the radial 
diffeomorphism constraint of the EYM theory. Thus there is practically no difference between 
the two theories concerning the number of degrees of freedom. In jHj a classification of all static, 
spherically symmetric solutions of the EYM theory with a regular origin was given; our aim is 
to prove a corresponding classification for the YMD theory. Naively one could expect that this 
should be an easier task for the YMD theory due to its simpler structure. However, life is not so 
simple and, although ultimately the result is the same, the proof for the YMD theory appears 
to be more difficult. The main reason is that in the EYM theory the 'area variable' r can have 
at most one maximum, its counterpart in YMD theory can oscillate. Apart from this subtlety 
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things in the EYM and YMD theory turn out to be largely the same. In both cases there are 
three different types of solutions with a regular origin. The generic one develops a singularity 
of the gravitational resp. dilaton field for a finite value of the corresponding autonomous radial 
variable. The second type is a countably infinite family of globally regular solutions differing by 
the number n of nodes of the YM potential W. Finally there is an oscillating limiting solution for 
n — * oo. Based on this classification an existence proof for the three different types of solutions 
for the EYM theory was given jHj, which can be easily adapted to the YMD theory. Futhermore 
the asymptotic scaling law of the parameters of the BK solutions for large n derived in [£] can 
be straightforwardly translated to the YMD theory. 

For the SU(2) Yang-Mills field W% we use the standard minimal spherically symmetric 
(purely 'magnetic') ansatz 

W*T a dx^ = W(R) {TxdO + T 2 sin Odip) + T 3 cos 9dip , (1) 

where T a denote the generators of SU(2) and R the radial coordinate. 
The action of YMD theory has the form 

Symd = / (\d^*P - l^ e ^F^F a ^)d A x , (2) 

where g and k are the gauge resp. dilaton coupling. Inserting the ansatz Eq. (J2) into the action 
and making use of the spherical symmetry we obtain the reduced YMD action 



S = dR [t^- + ^(W* + • (3) 




The dependence on g and k can be removed by the rescaling <j> — > <j)/n, R — > Rn/g and S — > Sgn. 
The resulting Euler-Lagrange equations are 

lBV)' = 2 e 2 *(iy' 2 + ( '^ 2)2 ). (4a) 
R z 

Using r = In R as a coordinate and introducing 

W 2 — 1 

r = Re-*, N = 1 - R6' , U = e^W T = , (5) 

r 

we obtain the autonomous first order system of Riccati type (the dot denoting a r derivative) 

r = rN , (6a) 

W = rU, (6b) 

N = 1 — N — 2U 2 — T 2 , (6c) 

U = WT + (N — 1)U , (6d) 

f = 2WU - NT (6e) 

supplemented by the constraint W 2 — rT = 1. These equations have a great similarity to the 
Eqs. (50) of [H] for the Einstein- YM theory. Clearly there is no analogue of the diffeomorphism 
constraint here. 

There is a kind of r-dependent 'energy' E = 2W 2 — (W 2 — l) 2 obeying 

E = 4(2iV - l)W 2 , (7) 
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which will be useful as a 'Lyapunov Function'. In addition we introduce some other useful 
auxiliary quantities e, / and g defined as 

e = £ = 2U 2 - T 2 , f = (1 - N) 2 + e and g = l-N-f, (8) 
obeying the equations 

e = A(2N -1)U 2 -2Ne , (9a) 

f = -2f + AU 2 , (9b) 

g = -g + (l-N) 2 . (9C) 

(9d) 

In general it does not seem possible to solve Eqs. (jBJ in closed form, yet there are some simple 
exeptions. One is the trivial vacuum solution W 2 = 1, (p =const.; besides there is an analogue 
of the extremal magnetically charged Reissner-Nordstr0m (RN) solution of the EYM theory 

W = N = l-~ 4> = -ln(l + ce~ T ) . (10) 
r 

2 Singular Points 

Before trying to explore the global behaviour of solutions with a regular origin, it is important 
to know the singular points of the system ©. Actually, there are two types of singular points, 
those attained for finite r and the fixed points for r — > ±oo. The first type of singularity occurs, 
if the r.h.s. of Eqs. © blows up at some finite value of r. As we shall prove later, the only such 
possibility is that ./V — * —oo and r —* 0. Due to the simple Riccati form of Eqs. © it is easy to 
find all their fixed points. There are first the f.p.s r = 0, oo and N = 1,W = ±1, U = T = 0. 
Furthermore there is the f.p. N = 0,r = T = 1,W = U = like for the EYM theory. All these 
f.p.s are of hyperbolic type and thus the application of the theory of dynamical systems [3 
provides theorems on local existence and the asymptotic behaviour near the singular points. 

Introducing suitable auxiliary variables it is also possible to treat the singularity with — > 
— oo (and r — > 0) as a fixed point. For that purpose it turns out to be convenient to use r as the 
independent variable and two auxiliary dependent variables 



obeying the equations 



K = r(l-N), A = WT+ (AT — 1)U (11) 



dr k — r 

T u = ~—> ( 12b ) 

dr k — r 

dr k — r 

d x _(L-M + V*-f)U + \ (12d) 

dr k — r 

P^l^Hl. , 12 e, 
dr k — r 



(12f) 



Applying Prop. 4 of [£j assuming k / we obtain finite limits Wo, Uq, ko, Ao, A*o and /o for 



the dependent variables. The finiteness of / implies kq = \ Wq — 1| and thus we have to require 
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Wq 7^ 1. From the behaviour N ~ 1/r it follows that r resp. r stay finite as r — > and hence 
—* oo. The original dependent variables behave as 

^ = W ° + 2(1 ^iV 2 ) ^ + ^ + ° (r4) ' (13a) 
U = ±W - (3|W 2 - 1|W 3 + ^zj)r + 0(r 3 )) , (13b) 

|W n 2 - II o 

iV = — !— 9 L + 1 + JVir + 0(r 2 ) . (13c) 



with arbitrary parameters Wq, W3 and N\. Actually r = is a regular point of Eqs. ()12|) as long 
as k > and thus this type of singular behaviour with N — ► — 00 is of a generic type. 

The case Wq = 1 requires special treatment. Unlike the EYM theory there are no solutions 
with a 'Schwarzschild type' singularity Wq = ±1, k ~ |Wq — l\/r+No with iVo 7^ 1 and solutions 
with Wo = ±1 automatically have a regular origin with U = and iV = 1. The f.p. treatment 
of this case proceeds like for the Bartnik-McKinnon solutions given in [Hj and we don't repeat 
it here. The linearization at the f.p. gives the solutions 

W = ±(l-be 2r + ce" T ) , (14a) 
N = l -46 2 e 2 * e 2r (14b) 
(/) = 4>o + 2b 2 e 2 ^°e 2T + de~ T , (14c) 

where b, c , 4>o and d are free parameters. For regular solutions, i.e. those running into the f.p. 
for r — > — 00 the coefficients c and <i have to vanish. Without restriction one may choose the + 
sign for W and <po = 0; then one gets W = 1 — br 2 + c/r + 0(r 4 ) etc. 

Similarly there is the fixed point W = ±1, £/ = and N = 1 with r — > 00 corresponding to 
asymptotically regular solutions. The asymptotic behaviour close to the f.p. is again dominated 
by the linear approximation 

W = ±(l + ce~ T + be 2r ) , (15a) 
N = l + de~ T , (15b) 
4> = 4>oo + de~ T . (15c) 

Regular solutions require now 6 = and normalizing to ifioo = one obtains W = ±(l + c/r + 
6r 2 + 0(r~ 2 )) etc. 

There is one more fixed point analogous to the 'RN fixed point' of the EYM theory given by 
W = U = N = 0, r = 1 and hence eft — * 00. For simplicity we call it again the RN fixed point. 
Putting r = r — 1 the linearization at this f.p. is 



with the solution 





f = N, 


(16a) 




N = —N + 2f , 


(16b) 




W = U , 


(16c) 




U=-W-U , 


(16d) 


w = 


C ie -isin(^r + 0) , 


(17a) 


u = 


de 2 sm(— r + — + 9) . 


(17b) 


f = 


C 2 e^ - C 3 ±e- 2r , 


(17c) 


N = 


C 2 e r + C 3 e- 2r . 
4 


(17d) 



The oscillating solutions running into the f.p. for r — > oo require C 2 = 0. For these solutions 
the term C^e~ 2r does not describe the actual asymptotic behaviour of f and N, because the 
true behaviour is determined by nonlinear terms involving W and U. Integrating the resulting 
inhomogeneous linear equations for f and N one finds (neglecting the non-leading term C^,e~ 2T ) 

r = C 2 e T + Cfe- T (--^-cos{V^T + 29) + ^| sin^r + 20)) , (18a) 

\ Zo 28 ' 



N ■ 



28 28 
'3 , /7T 2^/3 



: C 2 e T + Cfe" 7 ' (- cos(v / 3r + 20) + -y- sin(V3r + 29)) . (18b) 



In order to run into the f.p. for r — ► —00 it is necessary that Ci and C3 vanish. From the form 
of Eqs. © one sees that the vanishing of C\ implies W = 0. Since the 'exterior' (r > 1) part 
of the limiting solution, when the number of nodes goes to infinity, runs into the RN f.p. for 
r — ► — 00 it must correspond to W = 0. 



3 Classification 

The classification of solutions with a regular origin proceeds very similar to the one for the 
EYM system in [5]. We distinguish three different classes. The first one, which we call singular 
(Sing) are solutions becoming singular according to Eqs. (|13|) with r — ► for some finite value 
to. This class may be further subdivided into the classes Sing n (n = 0, . . .) containing singular 
solutions with Wq > 1 and W having n nodes and the class Sing^ for singular solutions with 
Wq < 1. The second class are the globally regular solutions (Reg) reaching the f.p. with r = 00 
described in Eqs. (|15|) . Again the class Reg can be subdivided according to the number of zeros 
of W. Finally there is a third class, the oscillating solutions (Osc), running into the RN f.p. 
with r = 1. Our aim is to prove that any solution of Eqs. © with a regular origin belongs to 
one of the classes Sing, Reg or Osc. This is made plausible by the observation that the dilaton 
(ft is a monotonically increasing function of r. Thus it can either stay bounded or diverge to 
+00; the latter can either happen for some finite To or for r — * 00. The difficulty is to show 
that these alternatives lead precisely to the singular points described in Eas. ()13|) . (|15|) and 1)17(1 
corresponding to the three classes introduced above. The essential part of this claim are the 
equivalences stated in Props. (1) to (3), whose proof will be based on a series of Lemmas (some 
of which can be found in the Appendix). In what follows we always consider solutions of Eqs. (|fi|) 
starting from a regular origin with W = 1. 
Lemma 1: 

i) The function W can have neither maxima if W > 1 or > W > — 1 nor minima if 
W < -1 or < W < 1. 

ii) The functions (ft = 1 — N , f and g are hounded below and if any of them is positive for some 
To then it stays positive for all t > tq. In particular they are non-negative for solutions 
with a regular origin. 

in) For finite r the functions U and W are finite as long as N is finite. 

Proof: i) and ii) are trivial consequences of Eqs. © and @ and the fact that 1 — N, f and g 
are positive close to the origin, if they are regular at r = 0. 

In order to prove iii) we remark that lnr and U obey linear equations and stay obviously 
finite as long as W and iV are finite. Suppose now that iV is bounded and W is unbounded, then 
because of property i) \W\ actually tends to 00. Following closely an argument put forward in 
|Hj Prop. 5 one finds that W diverges for some finite tq like W = ±v2/ (to — T ) + 0(1). Plugging 
this into Eq. ©) leads to a divergence of N contrary to our assumption. 
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Lemma 2: Suppose N — ► — oo at tq then W and K = r(l — N) have finite limits Wq and kq at 
to with Wq / ±1 and ko > (and as a consequence r(ro) = 0). 

Proof: First we want to show that W is bounded. Assume the contrary and without restriction 
W > 1 and U > 0. If we can show that U/W is bounded we get that also W /W = rU/W is 
bounded and thus InW is bounded. Putting z = W(l — N) — UT one gets 

z = r(l- N)U-z > -z , (19) 

This implies that z is bounded below and therefore since T —* +oo we may assume— z/T < 1 
and estimate 

dr\w) V 'W W 2 V W 2 WTWJ K J 

The r.h.s. is negative, if U/W > 2 and thus U/W is bounded implying the boundedness of W. 
Then also U and W are bounded implying the finiteness of W(tq). 

Next we want to show that k has a finite positive limit. Putting / = r 2 f we get from Eq. © 

'f = 2(N -l)f + 4W 2 . (21) 

Since J T ° Ndr < and W 2 is integrable we get from Lemma El that / is bounded and hence 
also k. The latter implies r — > and thus J T ° -/Vdr = — oo. Using once more Lemma IT71 we get 
/ —> 0. Since the boundedness of U = W/r yields W — > we get k 2 — >• (Wq — l) 2 . 
In order to show that k stays away from zero we use again g > to estimate 

k=(N-1)k + 2rU 2 + rT 2 >-k + ArU 2 > -k. (22) 

This inequality shows that hiK is bounded below and consequently Wq ^ ±1. 
Lemma 3: If N < —1 at some point t\ then N tends to — oo for some finite tq > T\. 
Proof: From g > we get T 2 > 2£/ 2 - N(l - N) and from Eq. (j^) 

iV(r) < 1 - N - 4U 2 + iV - iV 2 < 1 - N 2 < for r > n (23) 

and thus iV stays below —1. Suppose N(t) < —1 — e for r > t\ with e > then iV < — 2e. 
Hence iV is unbounded from below. Using Lemma IT51 we get — > — oo for some finite To- 

Lemma 4: IfW 2 > 1 and WU > at some point then the function N tends to —oo for some 
finite tq. 

Proof: The proof follows essentially the one given in [Sj . From 

f = 2WU-NT, U = WT+(N-1)U, (24) 

we get 

4- In \TU\ > 2y/2\W\ - 1 > 2V2 - 1 , (25) 

It follows that for any c > we can find some n such that 2y/2\TU\ > c for r > n and hence 
T 2 + 2C/ 2 > c. Eq. ((SJ3) then yields N < 1 — c — N implyig that N eventually becomes arbitrarily 
negative for large enough r and thus Lemma 01 may be applied. 

Lemma 5: r cannot have a maximum with r max > \/2 for W 2 < 1. 

Proof: At a maximum of r we have N = Q and iV < 0. On the other hand g > implies 
T 2 > 2U 2 for iV = and thus 

iV>l-2T 2 >l-il>0. (26) 
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Lemma 6: If In r is bounded for r — > oo then A?~ — > 0. 

Proof: Lemmas n and |21 imply \N\ < 1. Suppose lim sup A?" > 0, then there exist some e > 
and a sequence of points Tj — » oo with Nfa) > e. Since iV < 2 we get A/"(r) > e/2 for 
Tj — e/4 < r < Tj, implying J°° N 2 dr = oo. Similarly the assumption liminf Af < leads to the 
divergence of J°° N 2 dr. On the other hand we have < g < (1 — TV) 2 < 4 and thus 

poo poo 

oo > / (g- N)dr = 2 (N 2 + 2U 2 - N)dr . (27) 

Since J°° Ndr = lnr+ const, is bounded this implies the boundedness of J°° N 2 dr contradicting 
the assumptions lim sup N > and lim inf N < and hence lim N = 0. 

Lemma 7: r — » oo for r — » oo implies N — * 1 for r — > oo. 

Proof: According to Lemma |1] we may assume W 2 < 1 and thus T -> for r — > oo. Thus for 
any e > there is some r e such that \T\ < e for r > r e . From g > we get 

A(i _ at) = at _ i + 2[/ 2 + T 2 < 2T 2 - (1 - AT) 2 < 2e 2 - (1 - N) 2 . (28) 
dr 

Together with 1 — N < 1 implied by Lemma |S] this shows that < 1 — N < 2e for r > r e + j"? • 
Hence — > 1 for r — > oo. 

Proposition 8: The following are equivalent 

i) <p — > oo for some finite tq 

ii) r — > for some finite tq 
Hi) N — > — oo for some finite tq 
iv) The solution belongs to Sing 



Proof: Obviously iv) implies i) and i) implies ii). 

ii) hi): Since lnr diverges must be unbounded from below and thus iii) follows from 
Lemma El 

iii) => iv): We have to show that r — > and the functions W, U, k, A and /i have a finite limit 
at To. From Lemma|21we know that r — > 0, VF — > T^o 7^ ±1 and k — > ko > 0. Thus we have only 
to prove that the r.h.s. of Eqs. (|12b .d.e) are integrable. For that reason we have to show that 
U 2 and U 3 are integrable, which will be achieved using arguments put forward in 0. 

The boundedness of r and 1/k implies that |rA^| _1 and r e are bounded for any e > and 
consequently 

Nr e dr = -r e < 00 (29) 



c 



and hence 

1 



r 



— e dT < 00 . (30) 



Applying Lemma El we get that r e U obeying the linear equation 

(r e U)- = W ^_~ - + (AT - 1 + e)r e U (31) 

is bounded for e > 0. As a consequence \U n \ is integrable for any n > 0. According to Lemma lT71 
this implies that [i and consequently A and U have a limit. 

Proposition 9: The following are equivalent 
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i) <p — ► (j)^ < oo for t — > oo 

ii) r — > oo for t — ► oo 

iii) The solution belongs to Reg 

Proof: Obviously iii) implies i) and i) implies ii). 

ii) iii): From Lemma we know that TV" — > 1 for r — > oo. Thus J5 is asymptotically 
monotonously increasing. Suppose E — > oo, then \W\ — ► oo. Yet, this is not compatible with 
the boundedness of W . Hence W must be bounded and in fact tend to zero for r —* oo according 
to Lemma EH Since E has a limit also has a limit, which must be a f.p. of Eqs. ©• The f.p. 
with W = is, however, excluded, because E is asymptotically increasing and thus cannot tend 
to its infimum —1. Therefore the solution belongs to Reg. 

Proposition 10: The following are equivalent 

i) 4> °° f° r T 00 

ii) | In r \ < c < oo for all r 

iii) The solution belongs to Osc 

Proof: Obviously iii) implies i). 

i) => ii): From LemmaEland Prop.|Hlwe conclude that r must be bounded and thus J°° Ndr 
is bounded from above. Eq. (|27|) also shows that J°° Ndr is bounded from below and thus In r 
has also a positive lower bound. 

ii) =>- iii): Lemma |H1 implies iV — ► for r — > oo. Now we can use Lemma to conclude that 
E has a limit and VF — > 0. Thus also has a limit, which must be a f.p. of Eqs. ©• Eq. (j2t) 
is only compatible with iV — > 0, if the f.p. is = and r = 1. 

Putting together the Props. IH1 FH and fTTH with the possible behaviours of 4> discussed at the 
beginning of this section we obtain the 'Classification Theorem' 

Theorem 11: Any solution of Eqs. |^|) with a regular origin belongs to one of the three classes 
Sing, Reg or Osc. 

4 Topology of 'Moduli Space' and Existence Theorem 

The method to prove the existence of at least one globally regular solution for each number of 
nodes and a corresponding limiting solution with infinitely many nodes used in [Hj can be almost 
literally translated to the present case. 1 The proof is based on an analysis of the phase space as 
a function of the parameter b determining the solutions with a regular origin. While the generic 
singular solutions (i.e. solutions in Sing) correspond to open intervals of b space, the b values 
for regular solutions are isolated points accumulating at the value(s) for the limiting solution(s). 
As a first step we will study, what happens for very small and very large values of b. The result 
is the same as for the EYM system. 

4.1 Small b 

Proposition 12: Ifb^O is small enough the solution with W\ r= o = 1 belongs to the singular 
class Sing for b < or to Sing 1 for b > . 

The numerical analysis yields exactly one regular solution for each node number and correspondingly one 
single limiting solution 
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Remark: In view of Lemma ^ the restriction for b to be small is unnecessary for b < 0. 
Proof: The proof runs completely along the lines of the one in Rescaling r — > \b\~zr and 
U — > \b\~2U we obtain from Eq. (JHt) 

N = 1 - N - \b\ (2U 2 + T 2 ) (32) 

and the b independent boundary condition lim T ^_ 00 y = =p2. For 6 = the solution is N = 1 
and thus r = e T . As was shown the resulting solution W of the pure YM system diverges like 
W ~ ±\/2/( r — t) for some finite r resp. r_|- depending on the sign of b. The values of r± have 
been determined in jS] numerically as r + =~ 5.317 and f_ ~ 1.746. For small 6 we obtain a 
small perturbation of this solution as long as |6| Jl OQ {2U 2 + T 2 )dr' <C 1. This condition still 

holds for \W\ S> 1 and |6|~2> ^> 1, if b is small. Thus the Prop, follows from Lemma|IJ 

Before we proceed to the case of large values of b we shall derive the asymptotic behaviour of 
Wq as r runs back to zero, similar to the one given in |TJ. This proceeds in several steps. First 
we integrate the Eq. (JBfc) for N from N = 1 to N = using the pure YM solution becoming 
singular at r = f . We obtain some value r e < f from 

1= l/ Te dr' 8\b\ 

fl /-oo(f-T') 4 3f2(f-r e )3- 

This gives \W e \ ~ ; ~ 3 a | fe| ~ » /2. Here we have neglected the term 1 — N in Eq. (JHt) 

and pulled r out of the integral; this yields an error of 0(1) to W e , small by one order in \b\ 3. 
Correspondingly we get r e = f±. 

The next step is to integrate Eqs. @ from N = to N = — oo and to determine the value 
of Wq. Since rN is O(l) for N — > — cxd and r e = 0( | i>| — 2 ) we rescale r — > r\b\5 and similarly 
the dependent variables W — > Tyl6|~2 etc. Keeping only leading terms as 6 — > we obtain from 
Eqs. © 

r = rN , (34a) 
VF = rI7 , (34b) 
N = -2U 2 - T 2 , (34c) 
U = WT + NU , (34d) 
T = 2WC7 - iVT (34e) 

with the constraint rT — W 2 = 0. Due to the scaling the boundary conditions at r e are changed 
to r — > —00 and r = f±,W = U = T = N = 0. In order to perform the limit r -> we 
introduce the analogue of the variables used in Eq. ((12)1 putting k = —rN, A = WT + NU and 
/ = iV 2 + 2U 2 — T 2 . The Eqs. (|34[) imply / = and hence / = 0. Using the new variables we 
obtain 

r = —k , (35a) 
iy = W7, (35b) 
U = A , (35c) 
A = (3iy 2 -4U 2 )U (35d) 
K = 4rt/ 2 . (35e) 

These equations have to be integrated from the highly degenerate f.p. W = U = k = A = 0,r = 
f-t attained at r = —00. In order to lift the degeneracy we use a 'blow up' in the direction W 
12 introducing 

d = w^ R = w and ~ x = w (36) 
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and using In \ W\ as new independent variable. Thus we obtain the equations 

dr W 2 R , „ , 

w m—fo< (37a) 

W§f = ~-2U, (37b) 
dW rU 

„ r dX 3 - 3rA - 4iy 2 t7 2 

w w= ; (37c) 

w w = - {4 ° + 3S) • (37d) 

The f.p. has now moved to the point 

1 2 /2 1 

W = 0, r = f±, f/ = ±^^, k = = F ^ : — and A = — . (38) 

V2> ± 3r ± r ± 

Linearisation at this f.p. yields exclusively negative eigenvalues and thus the degeneracy has 
been removed. At the same time this shows that there are no adjustable parameters at the f.p. 
as to be expected. 

Numerical integration of Eqs. (|37j) from the f.p. to r = results in a finite value of W 
proportional to f±. Taking into account the scaling factor \b\~2 we find 

W ~ 0.89369| b\~ 3 f± . (39) 

Fig- CI shows a plot of numerically obtained data for Wq; the dashed lines represent the values 
of Eq. (|39|) for the two different signs of b. 

Next we turn to solutions with large values of b. Again the situation resembles very much 
the EYM case. 

4.2 Large b 

Proposition 13: If b S> is large enough the solution with W\ r =o = 1 belongs to the singular 
class Sing^. 

Proof: We put r = f/b and W = 1 + W/b. Keeping only leading terms we obtain from Eqs. (jfjj) 

f = Nr;, (40a) 

W = f U , (40b) 

f = 2U — NT , (40c) 

U = T + (N- 1)U , (40d) 

N = 1 - N - 2U 2 - T 2 (40e) 

with the constraint 2W — FT = and the b independent boundary condition lim-j-^-oo U- = =f 2. 
The combination z = N — 1 + TU obeys the simple equation z = —z. In view of the initial 
condition lirriT-^-oo z = we get z = 0. This allows us to remove from the T, U system to 
obtain 

f = U(2 + T 2 )-T, (41a) 
U = T{l-U 2 ). (41b) 

with the boundary condition liniT-^-oo U = lim r _ + _ 00 T = such that linv^-oo ^ = 1. It is 
straightforward to analyze the flow of this 2d system (Fig.|51shows a plot). The point U = T = 
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Figure 1: Numerical data for Wo(b); the dashed lines represent the values of Equ. I)39|) for the 
two signs of b 
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Figure 2: The phase space of Eqs. (|41|): on the dashed curve U = 2+t' 1 ^ e orm ts are vertical 
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is a hyperbolic f.p. with the eigenvalues —2 and 1. No orbits can cross the lines U = ±1. The 
solution with the relevant boundary condition from above corresponds to the separatrix for the 
eigenvalue 1. For 1 > \U\ > \T\/(2 + T 2 ) the orbits are monotonously approaching one of the 
lines U = ±1 for \T\ — > oo. This shows that eventually N < —1 even taking into account the 
correction terms of 0(l/b) in Eqs. ©• Applying Lemma |3] together with W = 1 — 0(1/6) proves 
the proposition. 

4.3 Existence Theorem 

As for the EYM theory one can characterize the neighbourhood of the sets Reg n and Osc. 
Proposition 14: Given b n E Reg n for any n then all b ^ b n sufficiently close to b n are either 
in Sing n or Sing n+1 . 

Proposition 15: Given boo £ Osc and some no then all b ^ boo sufficiently close to boo are 
either in Sing^ or in U„> no (Reg n U SingJ. 

Equipped with the knowledge, what happens for large and small values of b one proves 
Theorem 16: 

i) The sets Reg n and Osc are all nonempty, i.e., for each n = 0,1,2,... there exists a 
globally regular solution with n zeros of W for at least one b n £ Reg„ and there exists an 
oscillating solution with N > for all r and r — > 1 for r —* oo for at least one boo £ Osc. 

ii) The union [J Reg n has accumulation points that are contained in Osc, i.e., there exists 

n>0 

at least one sequence of globally regular solutions and one oscillating solution Woo such 
that W n (r) — > Woo(r) for r < 1 and W n (r) — * for r > 1 for n — » oo. 

The proofs of the Props. El an dEl an d of the Theorem can be literally taken from [S] and 
will not be repeated here. 

Remark: The existence of at least one regular solution for any n was already proven in jSJ. 
As mentioned above, the numerical results show that there is exactly one regular solution for 
any n and correspondingly only one limiting solution. As in the EYM case there is however no 
uniqueness proof available or in view. 

5 Scaling law for large n 

In jHj a remarkably well satisfied asymptotic scaling law for the parameters of the regular so- 
lutions with a large number n of nodes of W was formulated. The derivation was based on 
the observation that for solutions with many nodes three distinctive regions could be observed. 
In an inner region the solutions are well approximated by the limiting oscillating solution with 
infinitely many zeros. This region extends between the origin and r = 1. Furthermore there 
is an asymptotic region for r > r n S> 0, where the solutions are close to the flat solution con- 
necting the f.p.s with W = dbl, U = and W = U = 0. In the intermediate region extending 
between r = 1 and r = r n the functions W and U stay small and thus the equations for W 
and U can be linearized on the metric background given by the extremal Reissner-Nordstr0m 
solution. Boundary conditions for these linearized YM equations are obtained by matching with 
the solutions obtained in the inner and outer regions. The same type of scaling law can be 
obtained here through more or less identical reasoning. 

On region I defined by < r < 1 the solutions are approximated by the limiting solution 
running into the f.p. W = U = N = 0. The corresponding behaviour near r = 1 is given by 
Eq. (fTTj) neglecting the C3 term. The (n dependent) coefficient C2 has to be positive to allow 
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the solution to reach region II with r > 1. By a suitable shift in r we can always achieve C% = 1. 
The phase 9 is adjusted such that ^-t + 9 = mir at the m th zero of W. 

In region II the Eqs. (Jf^d) for W and U are linearized in the background of the solution 
Eq, ([TT)|) . Surprisingly the solution is identical to that of [0] 

r n (r) = l + C7 2in e r ;, (42a) 

Wn(r) = e"3 T sin(^-r + 9) + C 2 , n e^ sin(^-r + | + 0) . (42c) 

In region III, where we have r » 1 and 1 — N <C 1, we take the flat solution connecting the 
f.p.s with W = dbl, ?7 = and = £/ = 0. In the region where W is small this solution can be 
approximated by 

rm (t) =c n e r , (43a) 
Win(r) = ±£i,„e5* sin(^f + - + 0) , (43b) 



with the normalization Win — > ±(1 — e T ). Again the phase 6 is adjusted such that 

2 

obtain 



+ | + 6 = —rmr at the last but m th zero of W. Matching m, W\\ with rm, Wm we 



C 2 , n e T = c n e f , (44a) 
C 2 ,n^ T = Cie3 f , (44b) 

AT /q 

— r + 61 = — f + 9 + rm , (44c) 
where n is the total number of zeros of W. Eliminating r and f we obtain 

C 2 , n = (7 ie 7S ( ^- n7r) = C 2 , e- n 7S (45) 

and 

cn = Ctffi = c e n 73 . (46) 
Since the coefficient C 2 (6) has to vanish for the limiting solution, i.6. 6 — &oo > gGt 

C 2 (6) = -Qfi^ib - b^) + 0((6 - 6oo) 2 ) • (47) 
Numerical integration of the limiting solution and its variation with respect to b yields 

9 « 1.562209 and C 2 0.835060 • (6^ - 6) , (48) 
while numerical integration of the flat YM equations yields [H] 

« 0.339811 and C x 0.432478 . (49) 
From Eqs. ()45I46() we obtain 

b n = b OQ - 1.04894 • e~ na , and c n = 0.213530 • e na , (50) 

with a = 4- ps 1.81380 and e a « 6.13371. 

v3 
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Table 1: Parameters b and c of regular solutions;numerical results versus asymptotic formula 
Eq. (JSUJ 



n 


^num ^asy 


Cnum C aS y 


1 


0.2608301456037 0.20848171 


7.525748e-01 1.309730e + 00 


2 


0.3535180998051 0.35161335 


7.320406e + 00 8.033504e + 00 


3 


0.3750018038731 0.37494861 


4.852149e + 01 4.927516e + 01 


4 


0.3787544658699 0.37875304 


3.014792e + 02 3.022394e + 02 


5 


0.3793733291287 0.37937329 


1.853097e + 03 1.853848e + 03 


6 


0.3794744134274 0.37947441 


1.137028e + 04 1.137096e + 04 


7 


0.3794908985808 0.37949090 


6.974585e + 04 6.974616e + 04 


8 


0.3794935863472 0.37949359 


4.278045e + 05 4.278025e + 05 



Table 2: Masses of regular solutions; numerical results versus asymptotic formula Eq. (|51(1 



n 




M asy 


1 


0.80380777208 


0.79578582 


2 


0.96559851724 


0.96670624 


3 


0.99432009439 


0.99457200 


4 


0.99907210998 


0.99911505 


5 


0.99984867329 


0.99985572 


6 


0.999975327358 


0.99997648 


7 


0.9999959774969 


0.99999617 


8 


0.99999934419598 


0.999999374 



Table 3: Quotients of parameters of regular solutions 



n 


Ab n 


Ac n 


AM n 


1 


4.56821367 


9.72714739 


5.70301662 


2 


5.78233316 


6.62825122 


6.05669973 


3 


6.07360547 


6.21331291 


6.12131340 


4 


6.12385583 


6.14668276 


6.13170021 


5 


6.13210160 


6.13582559 


6.13338085 


6 


6.13345216 


6.13404859 


6.13365394 


7 


6.13370605 


6.13376280 


6.13369695 
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In [S] also the asymptotic formula for the mass M n = 1 — fC^n was derived. The same 
formula is supported by our numerical data in the present case, but in contrast to the EYM 
theory we were not able to find a simple derivation. Putting in numbers yields 

M n = l- 1.25259 • e~ na , (51) 

Tables ^ an d El contain a comparison of the numerically determined, parameters b n: c n and M n 
with the asymptotic values computed with the formulas from above. Table |3] displays the 
quotients Ab n = , fe °°T fe " , Ac n = and AM n = = 1 ~^ In of the numerical data. All of them 

approach rapidly the value e 71 "/^. 

Appendix: 

Lemma 17: Consider a solution y of the linear differential equation y = a + by in some interval 
tq < t < 7~i with \a\ integrable. If 

c(r',r) = £b(r")dr" , (52) 

is bounded from above for tq < r' < r < t\ then y is bounded; if c(r', r) has a limit as r —* t\ 
then y{r) has a limit; if c(r', r\j = — oo then y(ri) = 0. 

Proof: All properties are implied by the explicit form 

y (r) = y(r )e c(ro ' r) + f a{r')e< T ' ^ dr' . (53) 

Lemma 18: Suppose y obeys the inequality y < a + by — y 2 . If a is bounded from above and b 
is bounded for r > To, then y is either bounded for all r > tq or diverges to —oo for some Bnite 
n > t . 

Proof: Let A, B be positive constants such that a < A and \b\ < B. We can estimate y < 
for \y\ > C = \2A + 2B and therefore y is bounded from above. Furthermore y monotonically 
decreases and (1/y)' > 1/2 for y < — C, and thus y — > — oo for some finite r\. 

Lemma 19: Suppose there is some open invariant subset I of the phase space of the system 
Eq. such that 2N — 1 has a definite sign in I. Then for any trajectory (2N — l)r 2 U 2 vanishes 
on all its limit points for r — > oo in I. 

Proof: Eq. (JJJ) shows that the function E is monotonous along trajectories in / and hence 
serves as a 'Lyapunov Function'. According to Lemma 11.1 of jS] E = (2A^ — l)r 2 U 2 vanishes 
on the limit points of the solution. 

Corollary: Suppose N — > 0, r — > tq ^ for r — > oo and W and {/ stay bounded, then the 
solution tends to a f.p. of Eqs. JS|) with {7 = and = or l^ 2 = 1. 

Proof: From Lemma H^l we know that W — > for r — ► oo. Since £J is bounded it has a limit 
for r — > oo and thus also W has a limit, which must be a f.p. of Eq. (JJJ) otherwise VF would 
not tend to zero. 
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